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Введение. Рассматривается кинематическая модель группового преследования множества целей на плоскости. 
Преследователи при достижении целей используют метод, подобный параллельному сближению. В отличие от 
метода параллельного сближения векторы скоростей преследователей и целей направлены произвольно. В 
методе параллельного сближения мгновенные направления движений преследователя и цели пересекаются в 
точке, принадлежащей окружности Аполлония. В групповой модели преследования множества целей 
преследователи стараются придерживаться сети прогнозируемых траекторий. 

Материалы и методы. В модели поставлена задача достижения целей преследователями в назначенные 
моменты времени. Она решается методами многомерной начертательной геометрии при помощи эпюра 
Радищева. Прогнозируемая траектория является составной линией, которая при передвижении цели 
перемещается параллельно самой себе. На плоскости проекций «Радиус кривизны — значение скорости» 
выводится допустимый диапазон скоростей преследователя в виде линий уровня (это прямые, параллельные 
одной из плоскостей проекций). На плоскость проекций «Радиус кривизны — время достижения цели» 
выводятся образы линий уровня скоростей. Ведется поиск точек пересечения образов линий скоростей с линией 
уровня назначенного времени. По линиям связи значения точек пересечения опускаются на плоскость «Радиус 
кривизны — значение скорости». По полученным точкам строим аппроксимирующую кривую и ищем точку 
пересечения с линией назначенной скорости. В результате получаем значения радиуса окружности при 
прогнозируемой линии траектории движения преследователя. 

Результаты исследования. По результатам проведенных исследований созданы тестовые программы и 
изготовлены анимированные изображения в системе компьютерной математики. 

Обсуждение и заключения. Данный метод построения траекторий преследователей для достижения множества 
целей в заданные значения времени может быть востребован разработчиками автономных беспилотных 


летательных аппаратов. 
Ключевые слова: многофакторный анализ, эпюр Радищева, цель, преследователь, траектория, радиус кривизны. 


Для цитирования: Дубанов, А. А. Анализ скорости и кривизны траектории в задаче преследования множества 
целей / А. А. Дубанов // Адуапсеа Епошеегие Везеагсв. — 2021. — Т. 21, № 3. — С. 275-283. 


№ рз://401.0го/10.23947/2687-1653-2021-21-3-275-283 
Финансирование. Работа выполнена при финансовой поддержке инновационного гранта Бурятского 
государственного университета в 2021 году «Управление четырехзвенным манипулятором по сигналам, 


полученным с нейроинтерфейса». 
© Дубанов А. А., 2021 





Информатика, вычислительная техника и управление 


[о®) 
х 
л 


БИр://уез-допза.га 


1 


А4уапсеа Епатеетиз Кезеагсй 2021. Т. 21, № 3. С. 275-283. 155М 2687—1653 





Апа[у$1$ оЁ (Пе $рее4 апд сигуа@иге ог Фе галестогу ш Фе рго ет о{ ригзито 
а $её о? багое($ 


А. А. Рирапоу (5) 54 
Вап7тагоу Вигуай Зе Ошуегзиу (Лап-Оде, Визз1ап Еедегайоп) 


24 аапдаБапоу(@та!.га 


Ттго4диспоп. А Ктетайс то4е!| о{ огоир ригзий оРа зе оЁ {агое{$ оп а р|апе 1$ сопз14егеа. Ригзиег$ изе а 1есфтаце 
пп|аг $0 рагаПе]| арргоасВ тео4 © асшеуе гоа|5. ОиПКе Фе рагаПе! арргоасВ тео4, Ше зрее4 уесог$ оЁ ригзиегз 
апа {агое{5 аге Чиецще4 ато ига Пу. шп Фе рагаПе] арргоасй те#о4, е шзатапеоц$ Фтесйоп$ оГ тоуетеп оЁ Ше 
рагзиег ап4 Ве {агое{ ищегзес а а рошЕ Беопоше ю Фе сие оР АроПопиа$. ш Ше этоир то4е о рагзише шире 
соа[5, Ве ригзиегз {гу фо аВеге 0 а пебмогК оЁргефдсае фга]лестютез. 

МаептаЁ апа Мешйо4$. Те то4е| 5е{5 е 1азК оЁ асшеуше гоа[5 Бу ригзиег$ аё дез1епаже ро ш Ите. ТЬ1$ 
рго ет 15 зо]уе4 Бу Фе шефо4$ ог шиаЧитепзопа| дезсириуе сеотехгу изше фе Каа1зВсВеу Фаогат. ТБе 
ргедсе4 Чаесогу 1$ а сотрозйе Ппе Фа тоуез рагаПе] о изе!{ \Веп Фе {агое тоуе$. Оп Фе рго]есйоп р|апе 
“Вага$ оЁ сигуабваге — зрее4 уае”, Ве регио 51е зрее4 гапое оЁ Ше ригзиег 1$ 415р]ауе ш Фе Юг оЁ1еуе] Ппез 
(Фезе аге згалее Ппез рагаПе| ю опе оЁ Ше рго]есйоп р|апез). Птасез$ о зрее4 1еуе] Ппез аге 415р]ауе оп Фе 
рго]есйоп р1апе “Ка4та$ оЁ сигуабиге — ите фо геасЬ Фе соа[”. ТВе зеатсВ Юг ро оЁ пиегзесноп оЁ Фе зрее4 Ппе 
ппасез ап Фе арроп\е4 ите 1еуе| Ппе 15 Беше сопдицеа. А1опе Фе соттишсайоп Ппез, бе уаез оЁ Ше 
пуегзесйоп рош аге 1о\еге №ю Фе р|1апе “Каз оЁ сигуайте — зрееа уае”. Озше Фе оМашеЯ ро, \е 
соп5гас{ ап арргохитайпе сигуе апа 1ооК Юг фе пщегзесноп рошЕё м Фе Ппе оЁ Фе азз1епеа зрееа. Аз а гези®, \уе 
се{ уашез ое гад$ оЁ Ше сибе ай Фе ргефсеа Ппе оЁ Ше га]есюгу оР Фе ригзиег. 

Кезий$. Вазе4 оп Фе гези{$ оЁ Фе сопдисеа гезеагсВ, {е${ ргостатз Вауе Бееп стеме4, ап апиизже пипасез Вауе 
Бееп та4де ш Фе сотрщег тафетайс$ зу$ет. 

Оёбсизяюп апа Сопсизбоп5. Ты тефо4 о{ сопзасйие фга]лесюпез оф ригзиег$ ю асшеуе а уапееу оЁ соа[$ а а 
отуеп ите уашез сап Бе ш 4етапа Бу деуе]орегз о{ ащюпотоц$ иптаппе4 аепа1 уе 1] ез. 


Кеунога$: пои птепзлопа| апа|уз1$, Вад1зВсВеу р[о%, фагоеф, ригзиег, {га]лестогу, га аз оЁ сигуайте. 


Еог спайоп: А. А. Оибапоу. Апа[у$1$ оЁ Фе зрее4 ап4 сигуахге оЁ Ве ‘талесюгу ш Фе рго ет о{ ригзише а 5е{ оЁ 
{агое{ 5. Адуапсед Епошеетио Кезеагсй, 2021, уа1. 21, по. 3, рр. 275—283. 6рз://401.0ге/10.23947/2687-1653-2021- 
21-3-275-283 


Еипат® трогтайоп: Фе тезеагсь 15 аопе \иф Фе Йпапсла! заррое Нот Вигуа Зайе ЧУшуегзиу ш 2021 (Фе 


шпоуаНоп этап “Сопёго! оГа Юиг-ПиК таптри|аюог Базе4 оп $1е1па15 гесеуе4 Нот Ше пеигопщег асе”). 


Введение. Рассмотрим модель расчета траектории преследователя на плоскости, где в каждый момент 
времени от преследователя до цели строится прогнозируемая траектория, которой постарается придерживаться 
преследователь (рис. 1). 
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Рис. 1. Групповое преследование множества целей 
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Кривые [4 (5), [5($), [3(5) состоят из сегмента дуги окружности и прямолинейного отрезка. В нашей 
модели радиус окружностей и есть ограничение по кривизне прогнозируемых траекторий движения 
преследователей" —- 

Задачу исследования в данной статье можно описать так: преследователи, двигаясь по 
однопараметрической сети прогнозируемых траекторий, должны достигать своих целей в назначенные 
моменты времени, в том числе и одновременно. Для решения выбраны методы многомерной начертательной 
геометрии с использованием эпюра Радищева. Однопараметрическая сеть (рис. 2) состоит из конгруэнтных 
линий параллельного переноса. Каждая линия — аналог линии визирования (это прямая, соединяющая 


преследователя и цель). 
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Рис. 2. Однопараметрические сети прогнозируемых траекторий движения преследователей 


Для модели выбрана комбинация окружности и прямой. Есть множество вариантов обеспечения 
ограничений по кривизне (спираль Корню, кубическая парабола, сопрягающаяся с прямой, ит. д.). Для 
достижения целей преследователи используют метод, аналогичный методу параллельного сближения [1-4]. 
Однако в нашем случае скорость преследователей направлена произвольно, а в методе параллельного 
сближения линии скоростей преследователя и цели пересекаются в точке на окружности Аполлония. 

В тестовой программе, написанной по материалам статьи, объекты передвигаются по квадрату 
[-30:30] х [-30:30]. Расчет идет в метрах. Исследования проводились для скоростей 20 м/с. Начальный 
радиус окружностей при прогнозируемых траекториях равен 2 м. 

Если в момент начала преследования преследователь находился в точке Р; с вектором скорости И;, то 
центр окружности С; радиуса т; будет находиться в точке: 
ть 

у, 
И 
Затем из точки положения цели Т; строится касательная к окружности (С;,т;). Совокупность 


и» 


касательной и окружности будет базовой линией прогнозируемой траектории движения преследователя [;($). 
Отметим, что в уравнении базовой линии из однопараметрического множества прогнозируемых траекторий 
параметризация производится от длины дуги. 

При новом положении цели Т; линия [;(5) смещается, оставаясь параллельной самой себе (рис. 3). 


'Банников А. С. Нестационарная задача группового преследования // Лобачевские чтения : мат-лы 5-й молодеж. науч. школы-конф. // Тр. 
математического центра им. Н. И. Лобачевского. Казань : Изд-во Казанского математического общества, 2006. С. 26-28. 

` Изместьев И. В. Задача преследования маломаневренных объектов с терминальным множеством в форме кольца // Геометрические 
методы в теории управления и математической физике: дифференциальные уравнения, интегрируемость : мат-лы междунар. конф. Рязань : 
Изд-во Рязанского гос. ун-та им. С. А. Есенина, 2016. С. 17—18. 
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Рис. 3. Итерационный процесс расчета траектории преследователя 


Допустим, 1-Й преследователь в момент Е; находится в точке РЕ; ‚ имея при этом прогнозируемую 
траекторию движения [., (5), соединяющую с текущим положением цели Те. В этом случае следующей точкой 
траектории преследователя будет Р; 

Р 


]+1° 


+1 — ЭТО точка пересечения линии Г. ;1($), которая соответствует положению цели Па В 
следующий момент времени &;.1 и окружности с центром РЕ; и радиусом и "ДЕ, ДЕ = Е... — 6. Такова модель 


< 4, 5 
построения траектории преследователя ’.. 


Рассмотрим задачу группового преследования, когда группа преследователей догоняет группу целей. 
Будем считать, что каждый преследователь Р; стремится достичь своей цели Т;, хотя у некоторых 
преследователей цели могут совпадать (рис. 1, 2). 

Причем преследователь Р; достигает цели Г; за определенное время Ё;, двигаясь с определенной 
скоростью И;. Для одновременного достижения целей необходимо равенство всех &; определенному значению. 

Из рис.1 видно, что для изменения длины базовой линии можно менять радиус касательной 
окружности. Касательная вводится, чтобы преследователь мог плавно перейти на прямолинейную траекторию. 
Если бы это было так, то задача была бы сведена к преследованию методом параллельного сближения. 

Начальная скорость преследователя направлена произвольно, что позволяет задействовать метод 
параллельного сближения с соблюдением ограничений по кривизне (рис. 2). Для этого используется составная 
базовая линия, которая при движении цели перемещается, оставаясь параллельной сама себе, происходит 
плавный переход к методу параллельного сближения с соблюдением ограничений по кривизне (рис. 3). Рис. 2 
дополнен ссылкой на анимированное изображение, где можно посмотреть плавный переход к параллельному 
сближению. 

Цель данной статьи — описание метода, при котором преследователь достигает цели в назначенное 
время из допустимых значений. Можно рассматривать и одновременное достижение целей группой 
преследователей [5-9]. 

Материалы и методы. По результатам исследований разработана тестовая программа одновременного 
достижения целей преследователями, которую можно посмотреть на ресурсе автора. Предложенный алгоритм 
реализует итерационную схему расчета траектории преследователя (рис. 3). 

В модели принята зависимость для преследователя Р, который достигает цели Т за время &: 

Е = Е(Р., Те, пр, Пт, Ир, Ут, В). 
Здесь Р., ТГ‚ — координаты точек положения преследователя и цели в момент начала преследования; Пр, Пт — 
единичные векторы направления движения преследователя и цели в момент начала преследования; Ир, Ит — 
модули скоростей преследователя и цели в процессе преследования; А — радиус окружности, смысл которой 
показан на рис. 1, 3. 

Фактически в модели подсчитывается число шагов, за которые преследователь достигает цели. При 
известном дискретном промежутке времени число шагов можно сопоставить с реальным временем. 


*Дубанов А. А., Севээн А.-К. Э. Кинематическая модель метода параллельного сближения : св-во о гос. регистрации программы для ЭВМ / 
Бурятский гос. ун-т им. Д. Банзарова. ВО 2020665641. № 2020664886 ; заявл. 20.11.2020 ; опубл. 27.11.2020. 

° Дубанов А. А., Севээн А.-К. Э. Моделирование траектории преследователя на поверхности методом параллельного сближения : св-во о 
гос. регистрации программы для ЭВМ / Бурятский гос. ун-т им. Д. Банзарова. ВО 2020666553. № 2020664893 ; заявл. 20.11.2020 ; опубл. 
11.12.2020. 

° Дубанов А. Догнать одновременно. Плоскость 1. ОВГ: Врз:// ууу. уоцеафе.сот/аюВ?у=7УМНМ№СЬ\Уто (дата обращения: 22.05.2021). 
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Если цель движется прямолинейно и равномерно, то зависимость времени достижения цели в уже 
начавшемся итерационном процессе можно считать функцией от двух переменных — модуля скорости 
преследователя и радиуса кривизны окружности: 

= Е(УЬ, В). 

В модели считается, что преследователь движется с постоянной скоростью Ир, однако ничто не мешает 
нам менять значения модуля скорости, как и радиуса кривизны. Допустим, что модуль скорости принимает 
дискретные значения из ряда Ур‚,1 Е [1: М], а радиусы окружностей на рис. 1, 3 принимают значения №,, ] Е 
[1: М]. 

Для дальнейших исследований применяется эпюр Радищева, где используются координатные 
плоскости (В,И)и (В, Е) (рис. 4). 


Функции времени, за которое 
преследователь достигает цели, 
соответствующие скоростям 
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Рис. 4. Определение радиуса окружности на эпюре Радищева 


На рис. 4 представлено экспериментальное построение временных зависимостей &;; = ЕР (УР В;). 
Графики на плоскости (В, Е) показывают, как время достижения цели зависит от радиуса окружности К при 
фиксированном значении скорости Ур. 

В качестве одного из оптимизирующих факторов [10-11] на плоскости (В, Е) выбирается равенство 
Е = К, где &% — требуемое время достижения цели. Далее для решения нашей задачи на плоскости (К,И) в 
качестве второго оптимизирующего фактора выбирается равенство Ир = Ур,, где Ир, — постоянная скорость 
преследователя. 

В постановке задачи говорится о том, что модуль скорости преследователя является неизменным. Тем 
не менее построенный ряд значений скоростей необходим для расчета радиуса окружности Ку на плоскости 
проекций (К,И). 

По линиям связи на плоскости проекций (В,ИУ) находятся соответствующие точки пересечения с 
линиями уровня скоростей Ур, (рис. 4). По полученным точкам в тестовой программе выполняется 
полиномиальная регрессия и в итоге получаем функцию зависимости скорости преследователя от радиуса 
окружности, при которой происходит достижение цели за время Ц. 

Затем ищем точку пересечения функции Ур = }(К) с линией уровня Ир = Ир. Абсцисса точки 
пересечения Ау и есть искомый радиус окружности, при котором преследователь Р за время & со скоростью 
Ур, достигает цели Г. 

Расчет ведется при условии, что цель движется равномерно и прямолинейно. Если цель изменяет 
направление или скорость, то рассчитывается новый радиус окружности составной базовой линии (аналог 
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линии визирования метода параллельного сближения), ставится новое время достижения при прежней скорости 
преследователя. 

При равномерном и прямолинейном движении цели низший предел времени достижения фиксируется, 
когда скорость преследователя направлена в точку К на окружности Аполлония (рис. 5). Данное положение 
рассматривалось в работах Р. Айзекса [12], Л. С. Понтрягина [13], Л. А. Петросяна [14—16], Н. Н. Красовского и 
А. И. Субботина [17]. 





Рис. 5. Окружность Аполлония 


Окружностью Аполлония называется геометрическое место точек, отношение расстояний от которых 
до двух заданных точек — величина постоянная, не равная единице: |РК]|/|ТК]| = |Ур|/Ут| (фис. 5). 

При рассмотрении множественного преследования группы целей в тестовой программе делают 
предварительный расчет траекторий движения преследователей при заданных начальных параметрах. Из 
времен достижения целей для расчета одновременного достижения выбирается наибольшее время, и оно будет 
критерием для расчета траекторий остальных преследователей“. Этот момент проиллюстрирован на 
анимированном изображении, где три преследователя достигают двух целей одновременно”. 

На рис. 6 показано, как для одного из преследователей было установлено более короткое время 
достижения цели. Рис. 6 также дополнен ссылкой на анимированное изображение, где можно посмотреть 
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Рис. 6. Достижение целей в различное назначенное время 
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Результаты исследования. На рис. 7 приведены некоторые результаты многофакторного анализа в 
задаче одновременного достижения цели двумя преследователями. 
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Рис. 7. Результаты многофакторного анализа в задаче одновременного достижения цели двумя преследователями 





Цель движется прямолинейно и равномерно. Для каждого преследователя построили ряд допустимых 
скоростей. Допустимые значения радиуса окружности варьируются при помощи дискретной переменной АР 
(шкала АК на рис. 7). 

На плоскости проекций (АК, Е) строим однопараметрическую сеть линий. Каждая линия соответствует 
определенному значению скорости и выражает зависимость времени достижения цели от приращения радиуса 
окружности. На рис. 7 показана однопараметрическая сеть линий скоростей одного из преследователей. Для 
второго в тестовой программе многофакторного анализа построена аналогичная сеть. 

Для каждого преследователя выбирается первый оптимизирующий фактор [10-11], отвечающий за 
одновременное достижение: 

= 

Здесь & — наибольшее из времен достижения цели, если бы преследователи независимо догоняли цель 
при таких же начальных условиях. 

На плоскости (АВ,Е) ищем точки пересечения линий уровня Ё=Ц с линиями скоростей 
однопараметрической сети. Точки пересечения находятся при помощи встроенных процедур решения 
уравнений. В системе компьютерной математики Ма®САЛ это может быть процедура гоог{. Найденным точкам 
пересечения отвечают значения АВ и У на плоскости проекций (АВ,И). 

К полученным точкам на плоскости проекций применяется встроенная процедура полиномиальной 
регрессии и находится характеристическая кривая зависимости скорости от радиуса окружности составной 
базовой линии (рис. 1). 
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На плоскости проекций (АК,И) рис. 7 изображена такая же характеристическая линия зависимости 
скорости и для другого преследователя. Далее применяется второй оптимизирующий фактор И, = И. = И. В 
тестовой программе объекты движутся с одинаковыми скоростями. Встроенными средствами компьютерной 
математики находятся точки пересечения с линией уровня У = \. Этим точкам соответствуют значения АА: и 
АЁЕ>. 

При запуске итерационного процесса зафиксировано одновременное достижение цели двумя 
преследователями (аннотация к рис. 2). При этом были найдены значения приращений АК. и АР. к начальному 
радиусу окружности и заданы: 

— значение времени достижения цели (5, 

— модули скоростей движения И, с. 

В статье описан метод достижения группой преследователей множества целей с возможностью 
назначать время достижения. Одновременное достижение целей — это частный результат данного подхода, 
который развивает метод параллельного сближения. При реализации метода в пространстве следует добиться, 
чтобы векторы преследователя и цели находились в одной плоскости”. 

Рассмотрим случай преследования в трехмерном пространстве при следующих условиях: мы хотим 
свести задачу к методу параллельного сближения, но скорость преследователя направлена произвольно. В этом 
случае базовую линию прогнозируемых траекторий движения преследователя следует строить в плоскости, 
образованной линией визирования и скоростью преследователя. 

Следующий шаг преследователя — точка пресечения сферы с радиусом, равным шагу преследователя, 
и базовой линии, параллельно перенесенной так, чтобы один ее конец совмещался с точкой положения цели. 

Перейдем к вопросу нахождения окружности Аполлония и точки К в трехмерном пространстве. Сама 
окружность будет находиться в плоскости, образованной линией визирования и скоростью цели. Определим 
такие параметры окружности Аполлония, как центр окружности (точка 0), радиус окружности т, точка 
Аполлония (А) и точка К. Для этого учитывают вектор скорости цели, модуль скорости преследователя, 
положения преследователя и цели. Есть аналитическое решение этой задачи в плоской системе координат 
(рис. 5). Центр координат находится в точке положения цели. Вектор абсцисс будет единичным вектором вдоль 
линии визирования, соединяющей положения преследователя и цели. Вектор ординат будет перпендикулярным 
вектору абсцисс, но в плоскости, образованной линией визирования и вектором скорости цели. 

Обсуждение и заключение. Использованные в данной работе методы многомерной начертательной 
геометрии основаны на вариации модулей скоростей и радиусов кривизны окружностей. При этом по условиям 
задачи модули скоростей преследователей являются неизменными. 

В статье приняты во внимание результаты, достигнутые в работах [18, 19, 20]. 

Предложенный подход позволяет анализировать модули скоростей и направлений начального 
движения. 

Модель четырехмерного пространства, представленную в работах В. П. Болотова (гиперэпюр 
Болотова)", следует задействовать для анализа: 

— модулей скоростей, 

— радиусов окружностей, прилегающих к преследователям, 

— начальных направлений движения преследователей. 

Приведенные в статье результаты исследований могут быть востребованы разработчиками 
беспилотных летательных аппаратов, которые выполняют групповые согласованные задачи. Роль оператора 
наведения может быть сведена к указанию целей и контролю выполнения задач. 
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